Macroscopic quantum effects in a strongly driven nanomechanical resonator 
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We investigate the nonlinear response of a vibrating suspended nanomechanical beam on external 
periodic driving. The amplitude of the fundamental transverse mode behaves thereby like a weakly 
damped quantum particle in a driven anharmonic potential. Upon using a Born-Markovian master 
equation, we calculate the fundamental mode amplitude for varying driving frequencies. In the 
nonlinear regime, we observe resonances which are absent in the corresponding classical model. 
They are shown to be associated with resonant multi-phonon excitations. Furthermore, we identify 
resonant tunneling in a dynamically induced bistable effective potential. 
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c3 



-a 
c 
o 

(N 
> 

o 

o 
o 

-I— > 



i 

C 

o 
o 



X 



I. INTRODUCTION 

The ongoing progress in miniaturization of microscale 
devices allows nowadays to fabricate mechanical res- 
onators on the nanometer scalei. To realize them in form 
of transversely vibrating beams, lithographically pat- 
terned doubly clamped suspended beamsS^^iLSiSiiSiii 
are designed, but also suspended carbon nanotubes dis- 
play mechanical vibrations^ (see also Ref. Il3|) . Beyond 
applications as electrometers^, for detecting ultrasmall 
forces and displacement s] 5 ' 11 , or radio-frequency signal 
processing^, the nanomechanical devices also allow to 
investigate fundamental physical phenomena. In par- 
ticular, due to their small size, the crossover from the 
classical to the quantum regime is of interesli 1 ^* 1 ^. The 
quantum behavior will be due to a macroscopic number 
of particles whose coherence is disturbed by the inter- 
action with the environment. Promising progress in ap- 
proaching the quantum regime experimentally has been 
reported recentl}&2iii. 

In this work, we propose a straight-forward way to re- 
veal the quantum behavior of nanoresonators via the non- 
linear response of a nanobeam to an external ac-driving. 
Mechanical excitation frequencies of nanodevices have 
been measured via ac _ c j r iving^iSiSiL2i2iiS4iiii&. We con- 
sider the amplitude of the fundamental transverse vibrat- 
ing mode of a doubly clamped suspended nanobeam un- 
der longitudinal compression. The compression allows to 
control the degree of nonlinearity of the nanoresonator. 
Starting from a continuum mode l 14 ' 15 (see Ref. for the 
classical counterpart), the interacting field theory closed 
to the Euler buckling instability can be reduced to the 
dynamics of a quantum particle in an anharmonic po- 
tential. By exciting the beam to transverse vibrations, 
its non-linear response can be determined. The environ- 
ment is included by coupling the system to a harmonic 
Ohmic bath 1 ^. For weak driving, the response will be 
maximal at the eigenfrequency of the nanobeam and sim- 
ilar to that of the analogous classical system. However, 
as we show below, for strong driving, the non-linear re- 
sponse of the quantum system is qualitatively different 
from its classical counterpart which is the standard Duff- 



ing systemit This feature allows to separate both the 
classical and the quantum regime. In particular, we find 
distinct resonances in the dependence of the amplitude 
of the fundamental mode on the driving frequency. They 
can be interpreted as resonant multi-phonon excitations 
and are determined by avoided level crossings in the Flo- 
quet spectrum. In turn, we identify a separation of time 
scales and resonant tunneling in the dynamically induced 
bistable effective potential. 



II. EFFECTIVE SINGLE-PARTICLE MODEL 
STARTING FROM A CONTINUOUS MODEL 

We start from the normal mode description of an elas- 
tic rectangular beam of length I, width w and thickness 
d such that I ^> w > diiiSiii, which permits to con- 
sider only transverse displacements of the beam. A lon- 
gitudinal mechanical force controls the nonlinearity of 
the potential energy in the transverse direction. At low 
temperatures, the higher modes are frozen out and the 
fundamental mode can be treated independently. The 
resulting effective potential energy of the fundamental 
mode contains terms which are quadratic and quartic 
in the transverse amplitude x of the fundamental mode. 
The externally applied ac-driving can be included on the 
same footing. The driving strength can be tuned within 
the fundamental mode description to the regime where 
non-linear effects come into play but the higher modes are 
still negligible (see below). One ends up with an effec- 
tive single particle quantum mechanical time-dependent 
Hamiltonian for the fundamental mode amplitude x act- 



ing as a position operator, i.e. 
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+ -x 2 + -x 4 + xfcos(ujt) . (1) 



Here, m* is the effective mass of the beam. The param- 
eters in the potential are obtained as^ 4 *^ a = to*Wq = 



*— 2 

m lu, 



o y % ) an< ^ = , where the frequency of the 

fundamental mode is given by cJq = ir 2 d^/ Q / p / {I 2 \f\2) . 
Q is Young's elasticity modulus and p the mass density of 



and $ 



the material. The longitudinal force generates the strain 
e = (I — Iq)/Iq where Iq is the equilibrium length of the 
beam. At the critical value e c = — 7t 2 g? 2 /(12/ 2 ) (for a rect- 
angular beam), the system reaches a bifurcation point 
which is the well-known Euler instability. The effective 
potential for the fundamental mode then is purely quar- 
tic. Close to the Euler instability, both quadratic and 
quartic terms appear. We consider the case of a monos- 
table potential, i.e., a > 0, which implies that the strain 
remains below its critical values, i.e., e < e c . This sit- 
uation is easier to be realized experimentally compared 
to the bistable potential a < 0. The ac-driving occurs 
with amplitude / and frequency u>. An upper limit / ma x 
for the regime of validity of the fundamental mode de- 
scription is given by the first harmonic threshold, i.e., 
a^o/max < Sa>i ~ 3hu}(r*. We scale the Hamiltonian 
(|T|) with respect to the units in space and time, i.e., 
xq = \J hj (m*o;o) and to = 1/woj respectively, yielding 
the energy scale huo and the dimcnsionlcss nonlinearity 
parameter = e ^_ £ . 

We include the effect of the environment phenomeno- 
logically by a set of harmonic oscillators which are bilin- 
early coupled to the system with the coupling constants 
c^. The Hamiltonian is 

The bath is characterized by the spectral density 

/H = ^^^-^) = m* 7 ^e-^, (3) 
2 ' mjbjj 
j J J 

where we have chosen the standard Ohmic form with 
damping constant 7 and with cut-off frequency lo c . The 
total Hamiltonian is H(t) = Hs(t) + Hb- 

We focus on the case when the coupling to the bath 
is weak and use a standard Born-Markovian master 
equationi^iSS which reads 

P = -\ [H s (t ) , p] + C ie} (p) + £ noisc (p) . (4) 

The commutator describes the coherent dynamics while 
the bath acts via the two superoperators for relaxation 
and noise, respectively, 

£rel(p) = [x, \p,p]+] , 

£noisc(p) = -j:[ x AQ,P}] ■ ( 5 ) 
The operator 

y* 00 

Q= drK(T)x H (t-T,t) (6) 
Jo 

involves the integral kernel 

f /*°° f Tiuj \ 

K(t) = - duil(uj) coth ( — ) cos ujt (7) 

tv J \2k B T J 
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FIG. 1: Amplitude A of the expectation value (x(t)) for 
varying the driving frequency ui. Parameters are hsT = 
0.1fiw 0l/ 3 = 0.1,/ = 0.1fej /zo,7 = 0.005^0- Dashed line: 
Results of the classical Duffing oscillator at T = with the 
remaining parameters being the same. Inset: Amplitude A of 
the classical Duffing oscillator for varying driving frequencies. 
The driving strength / increases from bottom to top. 



with T being temperature and the position operator 
x n(t — t, t) in the Heisenberg picture. This master equa- 
tion can be solved numerically to obtain p(t) and the 
expectation value 

(x(t)) = tr(p(t)x) . (8) 

In the stationary limit, (x(t)} shows oscillations with the 
frequency being the external driving frequency u) (plus 
a phase shift) and with amplitude A. This amplitude 
is the quantity of interest that allows to compare to the 
classical counterpart. 



III. NON-LINEAR RESPONSE OF THE 
NANORESONATOR 

With the oscillation amplitude A at hand, we can in- 
vestigate the non-linear response of the nanoresonator 
on varying driving frequencies u). The classical analogue 
shows a well-known rich dynamics including chaos and 
driving induced bistability. Without noise, we obtain the 
Duffing oscillato»iZi2i. Its response on varying the driv- 
ing frequency can be calculated perturbativelyii and is 
shown in Fig. ^ (dashed line and inset). For weak driving, 
the standard Lorentzian resonance of a harmonic oscilla- 
tor occurs. For stronger driving, the system enters the 
nonlinear regime, the resonance curve bends over and an 
effective bistability is induced yielding a hysteretical re- 
sponse. If Gaussian white noise is added, the range of 
bistability and therefore hysteresis is reduced 2 ^ since the 
system escapes more easily from the metastable state by 
thermal hopping. 




FIG. 2: Amplitude A for varying the driving frequency u> 
and for different choices of (a) the nonlinearity coefficient f3, 
(b) driving intensity /, (c) temperature T and (d) damping 
strength 7. The other parameters are kept fixed according 
to P = 0.1, f = 0.1hco /x ,k B T = O.lfi^o and 7 = 0.005w , 
respectively. Dashed line in (a) and (d): same as in Fig.0 



For the quantum case, one can ask whether a signature 
of the induced bistability still appears in the response. 
For a truly bistable potential, tunneling leads to a reduc- 
tion of hysteresis 23 . Also, resonant tunneling through 
the potential barrier leads to characteristic steps in the 
hysteresis cycled. However, the potential in the Hamil- 
tonian Eq. Q is monostable. 

We have calculated the amplitude A of the steady-state 
oscillations for varying driving frequencies ui. The result 
is shown in Fig. ^ The characteristic profile consists of 
a series of peaks and dips. The resonances are sharper 
for higher frequencies. For lower frequencies, the broad 
peaks overlap strongly and lead to a shoulder-like profile 
which is similar to the classical result (dashed line in Fig. 

The locations of the resonances are dominated by the 
parameters of the undamped driven system, see Fig. [2^ 
and [2h- On the other hand, by changing the parame- 
ters characterizing the bath (temperature T and damp- 
ing strength 7) the shape of the resonances is modified. 
As can be seen in Fig.|2t and|2Ji, for increasing temper- 
ature and damping, the sharp resonances are smeared 
out and finally fade out when the quantum coherence is 
suppressed. This indicates that the resonances are re- 
lated to resonant multi-phonon excitations in the driven 
nonlinear system. 



IV. RESONANT MULTI-PHONON 
EXCITATION 

In order to show this, we exploit the periodicity in 
time of Hs(t) and calculate the quasienergy (Floquet) 
spectrum for varying driving frequencies ui 2 ^ 2 ^. To this 




CO / CO, 



FIG. 3: Average energies Ei (top), quasienergy levels e a ^ 
(middle) and amplitude A of the fundamental mode for vary- 
ing driving frequencies u. The remaining parameters are 
f3 = 0.1, / = OAhujo/xo and 7 = 0.005oj - 

end, we solve the equation 

[H s (t) - ih^]\Mt)) =eMt)), (9) 

where the quasienergies e a are defined up to a multi- 
ple integer of hw. This means that the state \<j^a\t)) = 
e mw *|(/) Q (t)) is also an eigenstate of the Floquet Hamilto- 
nian, but with the eigenvalue 

£a,n = + nhw . (10) 

Since the quasienergies do not allow for global ordering, 
we have also calculated the mean energies averaged over 
one driving period 

a>,n I Cot : n )> (11) 

n 

where the |c a ,„) are the Fourier components of the Flo- 
quet states 2 ^, which can be obtained as 

\c a ,„) = ^ J 2 J' U dte m ^ a (t)) . (12) 
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FIG. 4: Decay rate F for the slow dynamics of the amplitude 
A for approaching the steady state. Symbols: Solution from 
the numerical iteration of the Born Markovian master equa- 
tion, solid line: Smallest non-zero eigenvalue of the rate ma- 
trix of an improved Floquet Markovian master equation (see 
text). Inset: time-resolved dynamics of A for to = 1.167ljo 
(solid line). Fhe dashed line shows a fit to an exponential 
e~ Ft . Here, f3 = 0.1, f = 0.1hw /x , k B T = O.lfiwo and 
7 = 0.005^0- 



For the harmonic case (/3 = 0), the quasienergies are 
all degenerate for infinitesimal driving. The strong driv- 
ing and the nonlinear potential lift this degeneracy. As 
follows from the result shown in Fig. the quasiener- 
gies show avoided level crossings (and the mean energies 
show exact crossings) for particular values of the driv- 
ing frequency lu. To these (anti-)crossings correspond 
the particular resonances in the frequency-dependence of 
the amplitude A. The resonances occur when the (un- 
damped) quantum system absorbs an integer multiple of 
fvjj. The shape of the resonances is related to the splitting 
of the quasienergy levels at the avoided level crossings. 
For lower frequencies, the level splitting is quite large. 
Moreover, the avoided crossings are not well separated in 
this regime. This implies that the resonances are broad 
and strongly overlap yielding the shoulder-like behavior 
of the response profile. For larger oj, the avoided level 
crossings become well separated and the energy splittings 
decrease. Note that the resonances occur despite the fact 
that the thermal energy is larger than the quasienergy 
level splittings at the avoided level crossing. In turn, this 
leads to sharp resonances in the amplitude. The physics 
of avoided quasi-energy level crossings has been discussed 
in detail in Ref. 

V. DYNAMICALLY INDUCED BISTABILITY 
AND RESONANT TUNNELING 

The bistability of the steady state of the classical Duff- 
ing oscillator does not survive in the quantum system 
since it escapes the metastable state asymptotically via 
tunneling 24 . Nevertheless, we find signatures of bistabil- 



ity and tunneling if we consider how the steady state is 
reached. For this, we show in the inset of Fig. 0] the am- 
plitude A (local maxima of the vibrations) for increasing 
time (starting with the ground state of the undriven os- 
cillator as the initial state). We observe fast oscillations 
at short times. They decay on a time scale 7 -1 which re- 
flects intrawell relaxation in the metastable state. Then, 
starting from a metastable state at intermediate times, a 
slow exponential decay towards the asymptotically sta- 
ble state can be observed. By fitting to an exponential, 
we extract the decay rate T for various driving frequen- 
cies, see Fig. 0] crosses. The decay rate shows resonances 
at the same values of the frequencies where the avoided 
crossings of the quasienergy levels occur (see arrows). 
This slow dynamics can be identified as quantum tun- 
neling in a dynamically induced bistable effective poten- 
tial (see also Ref. l28j) . The peaks in r indicate resonant 
tunneling 24 from the meta- to the globally stable state 
both of which are dynamically induced. 

The rate for tunneling out of the metastable well can be 
obtained directly from an improved Floquet-Markovian 
master equation 2 ^. In the moderate rotating wave ap- 
proximation, the time-independent master equation 

Pap(t) = —Z (Eq - £p) Pap(t) + ^ C a 0,a'p>p a >i3'(t) 
n cifV 

(13) 

in the Floquet basis can be established. The transi- 
tion rates containing the influence of the dissipative bath 
reacts 

£aP,a'/3' — ^ (N aa >,n + Nppt^ n ) X aa ' >n Xp' p : - n 
n 

/3",n 

— Seta' ^ N a »pt tn Xpi a ii,— n X a »p, n . (14) 
a" ,n 

Here, the coefficients are given by 

7TL*^V£ 1 

N a/3 , n = N{e a ~e(3 + nhuj) , N(e) = ee/fcflT _ 1 

X a p. n = [ ' dte- inut ^ a {t)\x\cf> {t)). (15) 
27T Jo 

The time-independent rate coefficients (|14|l can be writ- 
ten as a matrix which can readily be diagonalized numeri- 
cally. We find a clear separation of time-scales where the 
smallest non-zero eigenvalue indicates a slow tunneling 
dynamics by which the stationary state is approached. 
The solid line of Fig.0]shows the smallest non-zero eigen- 
value of the corresponding rate matrix which was ob- 
tained by direct numerical diagonalization using a basis 
of 8 Floquet eigenstates. The result agrees well with 
that from the numerical iteration of the Born-Markovian 
master equation. Note again the analogy to resonant 
tunneling in a static double-well potential 2 ^. The role of 
the eigenenergies in the static case is now played by the 



5 



quasienergies e a determining the coherent dynamics, see 
Eq. (T3J. In both cases, the avoided (quasi-) energy level 
crossings are the origin of resonant tunneling from the 
metastable towards the globally stable state. Neverthe- 
less, the incoherent part of Eq. I|13fl is crucial to observe 
the resonant tunneling in this driving induced bistability. 

VI. CONCLUSIONS 

In conclusion, we have found resonances in the am- 
plitude of the fundamental mode of a nanomechanical 
resonator which is driven into its non-linear regime. 
They occur for particular values of the driving frequency 
and are explained in terms of quantum mechanical 
resonant multi-phonon excitations between quasienergy 
states. This response profile is generic for the quantum 
system and is absent in the analogous classical system. 
In addition, we have identified resonant tunneling in a 
dynamically induced bistable effective potential. Since 
oscillation amplitudes of nanomechanical resonators 
can be detected with currently available experimental 
techniques, this macroscopic quantum effect should be 
measurable in the near future. A promising approach 
seems to be to couple a superconducting single-electron 
transistor capacitively to the resonator— allowing an 
ultra-high sensitivity. Indeed, a fundamental frequency 
of 19 MHz has been reported^ implying that quan- 
tum effects occur below a temperature of a few mK. 



Most interestingly, the reported quality factors around 
Q = 40000 seem to be very promising since we predict 
quantum effects at even smaller Q = luq/j, see above. 
Finally, we note that our model of a driven anharmonic 
mono-stable quantum oscillator is a generic model which 
finds applications in various other fields of physics. One 
important application is related to the non-destructive 
readout for superconducting flux or charge qubitsSS. 
There, a dc-SQUID which is inductively coupled to 
the qubit is driven by an ac bias current such that 
the SQUID remains in its superconducting state (for 
this, the amplitude of the ac current has to be smaller 
than the critical current of the SQUID). The Josephson 
inductance of the SQUID carrying information on the 
qubit state is measured for varying the frequency of the 
ac-current. In terms of our model, the SQUID provides 
a sinusoidal potential for the superconducting phase 
acting as a single quantum mechanical particle. The 
particle is initially localized in one potential minimum. 
The ac bias current provides the external driving and 
for stronger driving, the particle experiences the non- 
linearity of the potential well. The distinct resonances 
in the nonlinear response which we have found might 
help to increase the efficiency of the read-out process 
since they are sharper than the common linear resonance. 
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